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Abstract 



£> , We consider the transition from a spatially uniform state to a steady, spatially- 

periodic pattern in a partial differential equation describing long-wavelength convec- 
tion [1]. This both extends existing work on the study of rolls, squares and hexagons 
and demonstrates how recent generic results for the stability of spatially-periodic 
patterns may be applied in practice. We find that squares, even if stable to roll 
perturbations, are often unstable when a wider class of perturbations is considered. 
We also find scenarios where transitions from hexagons to rectangles can occur. In 
some cases we find that, near onset, more exotic spatially-periodic planforms are 
preferred over the usual rolls, squares and hexagons. 



1 Introduction 



Pattern forming instabilities arise in a wide number of physical and chemi- 
cal problems. Model partial differential equations are used to try to capture 
the essential features of the observed transitions. In many interesting exam- 
ples such as Rayleigh-Benard convection and react ion- diffusion problems, the 
model equations are invariant under all translations, rotations and reflections 
in the plane and patterns arise at a transition from a trivial solution consist- 
ing of no pattern. Linear stability analysis of the trivial solution leads to a 
critical curve describing how the wavenumber for instability, k, depends on a 
parameter, fi, in the problem. For parameter values below the critical curve 
the trivial solution is stable. At a critical parameter, /i c , instability onsets at 
a critical wavenumber k c . 

Weakly nonlinear analysis is often used to try to predict the type of patterns 
observed once the trivial solution becomes unstable. Two aspects make this 
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particularly difficult: firstly, the rotational invariance of the problem means 
that instability to a single wavenumber gives instability to a whole circle of 
wavevectors. In other words, if the trivial solution is unstable to rolls then it 
is unstable to rolls with any orientation in the plane. Secondly, for \x > /i c , not 
just a single wavenumber but a whole band of wavenumbers is unstable. Of- 
ten this second problem is addressed by assuming that, sufficiently close to fi c , 
boundaries in any real problem will select out one particular wavenumber and 
modes with neighbouring wavenumbers will be suppressed. In the case of the 
first problem, a tacit assumption is often made that, for a given wavenumber, 
only a finite number of critical wavevectors are relevant. For example, four 
critical wavevectors oriented at | to each other (see figure 1(a)) are chosen or 
six critical wavevectors oriented at | to each other are chosen (see figure 1(b)). 
In both cases the critical wavevectors generate a periodic lattice of points, a 
square lattice in the first case and a hexagonal lattice in the second. Con- 
sequently, the circle of critical wavevectors is replaced by a finite set and a 
finite dimensional centre manifold exists for the problem, of dimension four 
in the case of squares and of dimension six in the case of hexagons. If critical 
wavevectors are used which do not generate a periodic lattice then there is no 
reason a priori why a finite dimensional centre manifold exists, since modes 
arbitrarily close to critical occur. While non-periodic cases have been consid- 
ered [2], their validity requires an additional assumption on the suppression 
of these near critical modes. 

Weakly nonlinear analysis using wavevectors on a square lattice or a hexagonal 
lattice as shown in figure 1, provide a framework for examining the relative 
stability of either squares and rolls or hexagons and rolls respectively. In both 
cases generic bifurcation equations have been derived using symmetry argu- 
ments [3-5]. A more complete stability analysis for rolls has been performed, 
for example by Brattkus and Davis [6], who consider the relative stability of 
two sets of rolls oriented at an arbitrary angle for a problem arising in crys- 
tal growth. Similarly, a more complete analysis can be performed for squares 
and hexagons by considering families of different square and hexagonal lat- 
tices (for specific examples see figure 2). This problem has a high degree of 
symmetry and using group theoretic arguments, Dionne and Golubitsky show 
that, for each lattice, additional branches other than hexagons, rolls or squares 
bifurcate as primary bifurcations [9,10]. In spite of the fact that the Fourier 
transform of the new patterns involves only one critical wavenumber, in phys- 
ical space they appear to have more than one lengthscale. Such "superlattice" 
patterns have recently been observed in the Faraday crispation experiment 
[7,8]. (For examples see figures 4 and 8 below.) 

Using symmetry arguments, Dionne et al. derive the generic bifurcation equa- 
tions for the families of square and hexagonal lattices and examine the stability 
of certain primary bifurcation branches in terms of the coefficients of the bifur- 
cation equations [11]. This stability analysis enables two types of statement to 
be made. Firstly, since each lattice problem corresponds to a subspace of the 



original unbounded problem, and since hexagons and squares each exist on a 
whole family of lattices, the stability of these planforms can be considered to 
a countably infinite number of perturbations. While this is not equivalent to 
completely determining the stability of squares and hexagons in an unbounded 
domain, it does considerably extend previous results. Secondly, each individ- 
ual lattice corresponds to either a square or hexagonal domain with periodic 
boundary conditions. For each lattice, the relative stability of the primary 
branches known to exist from [9] can be calculated. These results, contained 
in [11], have not as yet been applied to any specific partial differential equation 
and it is this issue we address here. 

In this paper, we re-examine the relative stability of spatially-periodic so- 
lutions to a partial differential equation considered by Knobloch [1]. This 
equation, 



f t = af- /iV 2 / - V 4 / + kV • | V/| 2 V/ 

+/3V-V 2 /V/- 7 V-/V/ + ( 5V 2 |V/| 2 , (i) 

describes a number of long- wavelength partial differential equations which 
arise in convection problems. For example, when k — l,/? = 5 = 7 = 0we 
recover the planform equation for convection in a layer between two poorly 
conducting boundaries [12], and when k = 1, 7 = 0, f3 — — ^-,<5 = — ^- equa- 
tion (1) models long- wavelength Marangoni convection [13]. Further examples 
are given in [1]. For equation (1) we demonstrate, that with relatively little 
additional analysis, we can derive all the coefficients necessary to apply the 
results from [11]. We thereby significantly extend Knobloch's stability results 
by inclusion of the additional perturbations. 

In section 2 we define the critical modes which generate square and hexagonal 
lattices used here and the resulting generic bifurcation equations. The deriva- 
tion of the coefficients of the bifurcation equations is given in section 3. Then 
in section 4 we discuss the results for two specific cases. In Case I we take 
7 = 0, k = +1. Knobloch called this Case B, the nature of our conclusions 
for his Case A are similar and we do not present them in detail. For Case I, 
provided (3^5, the coefficient of the quadratic term in the bifurcation equa- 
tions for the hexagonal lattices is nonzero. This quadratic term renders all of 
the primary solution branches for the hexagonal lattices unstable at bifurca- 
tion [14] and thus we restrict our attention to the square lattice bifurcation 
problems. We divide our discussion into two parts: in section 4.1.1 we consider 
the stability of squares and rolls in an unbounded domain by considering their 
stability on the whole family of square lattices; in section 4.1.2 we consider the 
particular example of long- wavelength Marangoni convection and show that 
different bifurcation scenarios can occur for different square lattices. In Case 
II we consider 7/A; 2 = 5 — (3, k = +1. This choice of parameters yields a degen- 
erate bifurcation problem for the hexagonal lattice since the coefficient of the 



quadratic term in the bifurcation equations is zero. Stable primary branches 
are therefore a possibility for all lattices and we consider both square and 
hexagonal types. We first discuss what can be deduced of the stability of rolls, 
hexagons and squares in an unbounded domain in section 4.2.1; then in sec- 
tion 4.2.2 we discuss the unfolding expected if the coefficient of the quadratic 
term is non-zero but sufficiently small. Finally, in 4.2.3 we discuss the spe- 
cific example of Marangoni convection for different hexagonal lattices. Our 
conclusions are summarised in section 5. 



2 Preliminaries 



In order to apply the analysis given in [11], we consider sets of eight or twelve 
critical modes whose wavevectors generate square or hexagonal lattices respec- 
tively, where the length of the critical wavevector is greater than the distance 
between neighbouring points on the lattice. A typical example is shown in 
figure 2(a) for the square lattice. This figure should be contrasted with figure 
1(a), which shows the wavevectors which are used conventionally in pattern 
selection studies. In both figures, a circle representing the critical wavevectors 
for the original unbounded problem has been superimposed on the lattice. In 
figure 1(a) this circle only intersects the lattice at four points and there are 
consequently four critical modes, whereas in figure 2(a) the circle intersects 
the lattice at eight points and hence there are eight critical modes. Sufficiently 
close to the critical value of the parameter, /i c , all other modes, represented 
by vectors not of length k c , will be damped. A family of finer and finer lattices 
can be constructed each with eight points on the critical circle. Each lattice 
can be encoded by a pair of integers, (m, n); for example, the lattice shown in 
figure 2(a) corresponds to the case (2, 1) i.e. Ki s is two squares of the lattice 
across and one up. The eight wavevectors consist of two sets of four wavevec- 
tors, (±Ki s , ±K 2 J and (±K 3s , ±K 4 J, that comprise squares and are rotated 
by an angle 9 S relative to each other. An alternative way to specify each lattice 
is therefore through the lattice angle 6 S , where, 

, / 2mn \ , . 

\m z + rr / 



and m > n > are relatively prime positive integers that are not both odd. 
Reducing the circle of critical wavevectors to four critical wavevectors, as 
shown in figure 1(a), is equivalent to changing the original unbounded domain 
to a box whose side length is -r- and applying periodic boundary conditions. 
Using eight critical wavevectors, illustrated in figure 2(a) for the case (m, n) = 
(2, 1), corresponds to changing the domain to a box of side length V ^ rn fc +TO and 
again applying periodic boundary conditions. 



In a similar way, a family of hexagonal lattices exists where the number of crit- 
ical wavevectors is twelve. These twelve wavevectors consist of two sets of six 
wavevectors, (±K lh , ±K 2h , ±K 3h ) and (±K 4h , ±K 5h , ±K 6h ), that comprise 
hexagons rotated at the angle 9h given by 

_i f m 2 + 2mn - 2n 2 \ .. 

Oh = cos — — ■ — 2T » (3) 

\ 2{m z — van + n A ) I 



where now m > n > y are relatively prime positive integers, and where m + n 
is not a multiple of 3. An example of the case (m, n) = (3, 2) is shown in figure 
2(b). For both the square and hexagonal lattices the requirement that m and 
n are positive integers ensures that the critical wavevectors generate a peri- 
odic lattice. This is necessary if the centre manifold theorem is to be invoked 
to formally justify the use of finite-dimensional bifurcation equations. The re- 
maining conditions on m and n ensure that each lattice angle corresponds to 
a genuinely different case [10]. 

For the square case, letting Zj be the complex amplitude of mode e*( Kj ' s ' r ) 
where r = (x, y) then the generic bifurcation equations take the form 

ii = Xzi + (ail^l 2 + a 2 \z 2 \ 2 + a 3 j^3| 2 + a 4 \z 4 \ 2 )zi + 0(|z| 5 ), 

z 2 = \z 2 + (a 2 \zi\ 2 + ai|z 2 | 2 + a 4 \z 3 \ 2 + a 3 |^ 4 | 2 )^i + 0(|z| 5 ), 

i 3 = Xz 3 + (o 3 |^i| 2 + a 4 \z 2 \ 2 + ai\z 3 \ 2 + a 2 \z 4 \ 2 )zi + 0(|z| 5 ), (4) 

z 4 = Xz 4 + (a 4 \zi\ 2 + a 3 \z 2 \ 2 + a 2 \z 3 \ 2 + a 1 \z 4 \ 2 )zi + 0(|z| 5 ). 

One recovers the bifurcation equations associated with the wavevectors given 
in figure 1 by restricting to the subspace z% = z 4 = 0. Equations (4) have six 
known types of primary branch which are listed in table 1 along with their 
stability assignments in terms of the coefficients. Note that rolls and squares 
are the same on all lattices, but the rhombs (rectangles), super squares and 
anti-squares take a different form depending on (m, n). For example, changing 
(m, n) changes the aspect ratio of the rhombs. An example of one of the super 
square solutions can be seen in figure 4. Further examples of the different 
planforms are given in [9-11]. 

For the hexagonal lattice problem, letting Zj be the complex amplitude of 
mode e* ( - Kj h' r - ) then the generic bifurcation equations take the form 



ii = Xz 1 + ez 2 z 3 

+ (&iN 2 + b 2 \z 2 \ 2 + b 2 \z 3 \ 2 + b 4 \z 4 \ 2 + b 5 \z 5 \ 2 + 6 6 N 2 )zi + 0(|z| 4 ), 
z 2 = Xz 2 + ez 3 z x 

+ {p2\zi\ 2 + h\z 2 \ 2 + b 2 \z 3 \ 2 + h\z 4 \ 2 + b 4 \z 5 \ 2 + h\z 6 \ 2 )z 2 + 0(|z| 4 ), 
z 3 = Xz 3 + eziz 2 



Table 1 

Signs of eigenvalues for primary bifurcation branches on the square lattice; a±, . . . , a 4 

are coefficients in the bifurcation equation (4). 



Planform 


Signs of non-zero eigenvalues 


Rolls (R) 

7 = ^(1,0,0,0) 


sgn(ai), sgn(a 2 - a ± ), sgn(a 3 - ai), sgn(a 4 - ai) 


Simple Squares (S) 
z = A s (l, 1,0,0) 


sgn(ai + a 2 ), sgn(ai - a 2 ), sgn(a 3 + a 4 - a 4 - a 2 ) 


Rhombs (Rh s i )m , in ) 
z = A Rh (0 s )(l, 0,1,0) 


sgn(ai + a 3 ), sgn(ai - a 3 ), sgn(a 2 + a 4 - ai - a 3 ) 


Rhombs (Rh s2 ,m,n) 

z = ^(0 S + §)(1, 0,0,1) 


sgn(ai + a 4 ), sgn(ai-a 4 ), sgn(a 2 + a 3 - a 4 - a 4 ) 


Super Squares (SS min ) 
z = A ss (l, 1,1,1) 


sgn(a 4 + a 2 + a 3 + a 4 ), sgn(a 4 + a 2 - a 3 - a 4 ) 
sgn(a 4 - a 2 + a 3 - a 4 ), sgn(ai - a 2 - a 3 + a 4 ) 
sgn(^o), where ^0 = 0(A 2 <i™ + } ) 


Anti-Squares (AS mi „) 
z = ^ 5 (1, 1,-1,-1) 


same as super squares, except no — ► — /io 



+ (& 2 N 2 + &2N 2 + feiksl 2 + h\z A \ 2 + 6 6 |,2 5 | 2 + h\z 6 \ 2 )z 3 + 0(|z| 4 ), 
i 4 = Xzi + ez 6 z 5 (5) 

+ (6 4 |zi| 2 + b 5 \z 2 \ 2 + & 6 |-2 3 | 2 + fexl^l 2 + b 2 \z 5 \ 2 + & 2 1 ^e 1 2 ) ^4 + 0(|z| 4 ), 

+ {h\zi\ 2 + b±\z 2 \ 2 + b 6 \z 3 \ 2 + b 2 \z 4 \ 2 + h^ 2 + b 2 \z 6 \ 2 )z 5 + 0(|z| 4 ), 

Zg = AZq + 6Z5Z4 

+ (b 6 \ Zl \ 2 + b 5 \z 2 \ 2 + b 4 \z 3 \ 2 + 6 2 |,2 4 | 2 + b 2 \z 5 \ 2 + hlz^ZQ + 0(|z| 4 ). 

The standard hexagonal bifurcation problem is recovered by restricting to 
the subspace z 4 = z^ = z Q = 0. Primary branches for equations (5) are 
listed in table 2. Note that rolls and hexagons are the same on all lattices, 
but the rhombs and the super hexagons take a different form depending on 
(m,n). Examples of one of the rhombs and one of the super hexagon states 
are given in figure 6(b) and figure 8 respectively. As for the square lattice, 
further examples of the different planforms may be seen in [9-11]. 



3 Calculation of the Coefficients 



The form of the bifurcation equations (4) and (5) is determined by the symme- 
try of the problem. However, the coefficients, and therefore the stability of the 



Table 2 

Branching equations and signs of eigenvalues for primary bifurcation branches on 

the hexagonal lattice; e, b±, . . . , b% are coefficients in the bifurcation equation (5). 



Planform and branching equation 


Signs of non-zero eigenvalues 


Rolls (R) 


sgn(fci), sgn(eA R + (b 2 - &i)A£), 




z = ^(1,0,0,0,0,0) 


sgn(-eA R + (6 2 - 61)^), 




= XA R + hA 3 R + 0{A R ) 


sgn(6 4 - 61), sgn(6 5 - 61), sgn(6 6 -61). 




Simple Hexagons (H ) 


sgn{eA H + 2{b 1 + 2b 2 )A 2 H ), 




z = 4 H (l, 1,1, 0,0,0) 


sgn(-eA H + (61 - b 2 )A 2 H ), 




= XA H + eA 2 H 


sgn(-eA H + (64 + &5 + b 6 - h - 2b 2 )A 2 H ), 




+ (b 1 + 2b 2 )Al + 0(AJ I ) 


sgn(-eA H + 0(A 3 H )). 




Rhombs (Rh hl>m>n ) 


sgn(6i + 5 4 ), sgn(6i-6 4 ) 




z = A Rh (e h )(l, 0,0, 1,0,0) 


sgn(/xi), sgn(/x 2 ), where, 




= XA Rh + (b 1 + b 4 )A Rh + O(A Rh ) 


Mi + M2 = (-26i - 26 4 + 26 2 + 65 + &6)^Bft, 






/ii/x 2 = -e 2 ^4| fe + (61 + 64 - &2 - 65) (&1 + &4 " 


-b2~h)A Rh . 


Rhombs (Rh h2>m>n ) 


sgn(6i + 5 5 ), sgn(6i-6 5 ) 




z = A Rh (9 h + f)(l, 0,0, 0,1,0) 


sgn(/xi), sgn(/x 2 ), where, 




= \A Rh + (b 1 + b 5 )A Rh + O(A Rh ) 


A*i + A*2 = (-264 - 26 5 + 26 2 + 64 + h)A 2 Rh , 






ixnx 2 = -e 2 A 2 Rh + (61 + 65 - b 2 - 64) (61 + 65 - 


~b2-h)A Rh . 


Rhombs (Rh fe3)m , in ) 


sgn(6i + & 6 ), sgn(6i-6 6 ) 




z = A Rh (e h -f)(l, 0,0, 0,0,1) 


sgn(/xi), sgn(/x 2 ), where, 




= XA Rh + (b 1 + h)A Rh + O(A Rh ) 


A*i + M2 = (-264 - 26 6 + 26 2 + 64 + h)A 2 Rh , 






VWi = ~e 2 A 2 Rh + (64 + 6 6 - b 2 - 64) (61 + 6 6 - 


-b2-h)A Rh , 




sgn(eA SH + 2(64 + 2b 2 + b 4 + 65 + & 6 )4s# ) 




Super Hexagons (SH^ n )** 


sgn(eA S n + 2(6i + 26 2 - 64 - 65 - b 6 )A 2 SH ) 




z = A SH (l, 1,1, 1,1,1) 


sgn(-eA S H + 0(A| H )), sgn(-eA 5if + 0(A 3 SH ))* 


= XA SH + eA 2 SH + (6i + 26 2 )^| H 


sgn( / t/i), sgn(^ 2 ), where, 




+(6 4 + 65 + h)A 3 SH + O(A^) 


Hi + Hi = -4:eA SH + 4(6i - b 2 )A 2 SH , 
IHH2 = 4(eA SH ~ (h - b 2 )A 2 SH ) 2 






-2((6 4 - & 5 ) 2 + (64 - foe) 2 + (65 - 


h?))A% H 




sgn(/u ), where, ^0 = 0{A 2 ^~ 1 '). 





* These two eigenvalues differ at 0{A 3 SH ). 

** Super triangles [8] have the same eigenvalues except sgn(/Uo) — » -sgn( / uo)- 



different planforms, depend upon the particular application. One important 
result in [11] is that, although high order terms in the bifurcation equations 
are required to find the relative stability of some planforms, much is fixed 
by the cubic order truncation. Therefore, in this section we determine the 
coefficients of the quadratic and cubic terms in the bifurcation equations (4) 
and (5) from the long- wavelength equation (1). Our approach is to infer these 
quantities from the branching equations for rolls, simple squares, rhombs, and 
simple hexagons, which are computed using perturbation theory as follows. 
We let 



f = ef + e 2 f 1 + e 3 f 2 + ... 
A* — Vc + e/-*i + e/i 2 + . . . , 

where \x is the bifurcation parameter and instability onsets at the critical value 
He = 2(— a) 1 / 2 , at the critical wavenumber k c = (— a) 1//2 . In turn, we take 



Rolls 

Squares 

Hexagons 

Rhombs 



c.c. 



fo = A R e ikx 

f = A s (e ikx + e lky ) + --- + c.c. 

/„ = A H (e ikx + e ik( - x+V3 y)/2 + e ik{ - x ~^ y)/2 ) + ■■■ + c.c. 

k 



(6) 



A Rh (e ikx + e ik{cx+sy) ) + • • • + c.c, where c = cos(0), s = sin(0) 



At 0(e 3 ) the solvability condition gives an equation for the amplitudes A R , As, Ah, A R y l 
respectively. That is, 



A R = \A R + C R A R + 0(A 5 R ), 
A s = \A s + C s Al + 0(A s ), 
A H = XA H + eA% + C H A H + 0(A%), 
A Rh = \A Rh + C Rh A Rh + 0(A Rh ), 



(7) 



where A = k 2 (fi — fi c ) and 



^^{8 + ^-S 

c «=-*K^ + '-isO(' + isO 



C H = C R -^U K+ h + 2 S -±)(2 S + ± 



C Rh (6) =C S - 4A; 4 cos 2 OIk + PS- 



Pe 



fc 4 (l-4cos 2 



8e 2 (l-2cos 2 #) 
fc 8 (l-4cos 2 #) 2 



Note that, as expected, Crh{\) = Cs- Also note that C^h diverges as 9 — ► | 
since when 9 — | there is resonance between e* fc:c and e * fc ( c:r + s J/). 

The cubic coefficients in the equivariant bifurcation equations (4) and (5) 
are readily expressed in terms of the branching coefficients Cr, Cs, Cr and 
Cru- For instance, if (4) is restricted to the simple squares subspace z = 
(^,^,0,0), we find 

A s = XA S + (a x + a 2 )A%. (9) 



Comparing equation (9) with the appropriate branching equation from equa- 
tions (7) gives Oi + a 2 = Cs- Similarly, by restricting to subspaces for rolls 
and rhombs for the square lattice problems (4), we find 



a i — Cr, a 2 — Cs — Cr, 

a 3 = C Rh (9 s )-CR, a 4 = C Rh (9 s + ^)-C R , (10) 

where 9 S G (0, |) takes on one of the discrete set of values (2). Similarly, re- 
striction to subspaces for rolls, hexagons and rhombs for the hexagonal lattice 
bifurcation problems (5) leads to 



i>i — Cr, 62 — ~(Cff — Cr), 64 — CRh(9 h ) — Cr, 

h = C Rh \9 h + —J — C R , b 6 = C Rh \9 h — —J — C R , (11) 

where 9h G (0, |) takes on one of the values in the discrete set (3). 

Note that the expressions for Cr, Cs, e, and Ch are given in [1] and, in that 
paper, are used to calculate the coefficients a\,a 2 ,b\ and b 2 - The remaining 
expression for Cru is the only additional calculation required to enable all the 
remaining coefficients in the bifurcation equations (4) and (5) to be found. 



4 Stability Results 



In this section we use the bifurcation equations (4) and (5) to determine the 
relative stability of the steady planforms which are given in section 2. The 
results depend on the parameters k, (3, 5, and 7 in the long-wave equation 
(1). They also depend on the size of the periodic domain through the lattice 
parameters (m, n). We restrict our discussion to two cases and, where possible, 
compare and contrast our results with those given in [1]. 



The task of establishing the stability of each planform is considerably eased by 
the fact that the eigenvalues determining the relative stability of each primary 
branch are derived in terms of the coefficients of the Taylor expansion of the 
bifurcation equations in [11]. The signs of the eigenvalues are summarized in 
tables 1 and 2 for the square and hexagonal lattices, respectively. The sign of 
the first quantity listed for each planform gives the branching direction; if this 
eigenvalue is negative (positive), then the branch is supercritical (subcritical). 
If e ^ then simple hexagons and super hexagons bifurcate transcritically; all 
other patterns arise through pitchfork bifurcations. We distinguish between 
the two branches of hexagons, denoted H + and H~, which satisfy Ah > 
and Ah < 0, respectively. Similarly, there are two distinct branches of super 
hexagons, denoted by SH^j n . Omitted from the tables are the zero eigenvalues 
associated with translations of the patterns, and also information about the 
multiplicities of each eigenvalue; this information can be found in [11]. Note 
that certain eigenvalues in tables 1 and 2 are not determined at cubic order in 
the bifurcation equations. For instance, the relative stability of super squares 
and anti-squares depends on a resonant term of 0(\z\ 2m+2n ~ l ). However, in 
this case, if super squares and anti-squares are neutrally stable at cubic order, 
then, generically, exactly one of the two states is stable. There is an analogous 
stability result for super hexagons and triangles [8]. 

In the case of the square lattice, the eigenvalues that depend only on a± and 02 
can be determined by considering the restricted bifurcation problem z% = z^ = 
in equation (4). Similarly, those results for the hexagonal lattice that depend 
only on e, b\ and 62 can be obtained by considering the simpler hexagonal 
bifurcation problem. In general, the signs of all the remaining eigenvalues are 
dependent on the choice of lattice. However, in the special case when e = 0, 
we find that those eigenvalues which are unchanged on permutation of 03 and 
04 are independent of 8 S , and those which are unchanged on permutation of 
64,65 and 6 6 are independent of 6h- This is due to the particularly simple 6- 
dependence of CjihiO) in equation (8) when e = 0. Also, in the case of the 
square lattices, results for angles close to | and | must be interpreted with 
care because of the singularity in Cnh{d) at — f that occurs due to resonant 
interactions. 

In each of the cases we discuss below, we evaluate the signs of the eigenvalues 
for each planfrom and determine if and where they change sign. We present 
the results in the form of bifurcation sets separating different regions of sta- 
bility and instability for the relevant patterns. Note that along the stability 
boundaries themselves, the bifurcation problem is degenerate and the bifurca- 
tion equations (4) and (5) are insufficient to locally determine the bifurcation 
structure. These degenerate points have been analysed by Knobloch [1]. 
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4-1 Case I: 7 = 0, K — +1 



In case I, solutions on the hexagonal lattices are unstable at onset unless /3 = 5 
i.e. e = in equations (5). Thus here we focus on the square lattice problems 
and defer discussion of the hexagonal cases to section 4.2. 



4-1.1 The stability of squares and rolls 

If we can show that squares or rolls are unstable on any one lattice then 
they must be unstable in the original unbounded problem sufficiently close to 
onset. From the hexagonal lattice result, we can immediately infer that rolls 
are unstable to hexagonal perturbations for (3^5. Here we show that squares 
are also unstable unless (3 is sufficiently close to 5. The bifurcation sets for 
rolls and for squares in the (/3, <5)-plane are presented in figures 3(a) and (b) 
respectively. 

The sign of the branching eigenvalue for rolls, sgn(oi), is always negative, 
so rolls always bifurcate supercritically. The sign of a 2 — a\ determines the 
relative stability of squares and rolls and the corresponding bifurcation line, 
ai = a 2 in figure 3(a), is identical to the line o(0) =0 given in figure 2(b) of 
[1]. The relative stability of rolls and the two rhombic patterns is determined 
by the signs of 03 — Oi and 04 — 01. Since 03 and 04 are dependent on 9 S , the 
precise position of the corresponding bifurcation curves a\ = a 3 and 01 = a 4 
depends on the choice of lattice: those shown in figure 3(a) are for the case 
(to, n) = (2, 1). Qualitatively the picture is the same in all cases corresponding 
to < 6 S < |. For 6 S > |, the picture is similar except a 3 — a± and a 4 — a\ 
switch roles. The region of stability is indicated by the shaded wedges in the 
{j3, 7)-plane between a\ = a 2 and 01 = 04 (01 = a 2 and a\ = 03 if 6 S > j). 
As the lattice angle, 9 S , approaches f or | the region of stability of the rolls 
is reduced to a narrower and narrower region occurring only for large \(3\ and 
\5\. At precisely — | (or the complementary |) the hexagonal lattice must 
be considered. 

In figure 3(b) we show the analogous bifurcation set for squares showing where 
each of the three expressions given in table 1 for the eigenvalues for squares 
change sign. The lines 01 = a 2 and a\ = — a 2 correspond to the lines a(0) = 
and Pat(0) = given by Knobloch. In his study he found the squares were 
preferred to rolls for the region between these two curves. However, this region 
is significantly reduced when instability to super square (or anti-square) states 
is included through the eigenvalue a 3 + a 4 — a x — a 2 . The position of the 
corresponding bifurcation line given by 03 + 04 = 01 + a 2 is again dependent 
on the value of the lattice angle. As the lattice angle approaches | or | the 
region of stability of the squares is reduced to a narrower and narrower region. 
Interestingly this narrow region always includes the line /3 — S, for which the 
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hexagonal problem is degenerate. Recall that, for hexagonal lattices, it is only 
in this degenerate case that stable planforms can exist at onset. Thus, in 
summary, we find that there are only stable squares or rolls when /3 ~ S, that 
is e a: 0. 

Knobloch also considered the case where (3 = 5, 7 7^ which he refers to as case 
A. In this case we find that the region of stability for the rolls and for squares 
is diminished to a narrow region about 7 = as the lattice angle approaches 
I or |. Since 7 = corresponds to the degenerate hexagonal problem, again 
we find that squares and rolls are unstable unless e ps 0. 



4-1.2 Stability of other planforms: the example of Marangoni convection 

As discussed in section 2, for each lattice, there are in fact six primary branches 
known to exist, any one of which could in principle be stable. Since many of 
the eigenvalues are dependent on the lattice angle 9 S , the precise region of 
stability for each state is dependent on the choice of lattice, i.e. on the size of 
the domain for a box with periodic boundary conditions. Consequently, for a 
given physical problem, which planforms are stable at onset, can be dependent 
on the size of the box. We illustrate this with the example of Marangoni 
convection which corresponds to (3 — — ^-,S = — ^- and 7 = 0. This lies 
within the region where squares are preferred to rolls in Knobloch's analysis. 
In contrast, we find that one of the following scenarios occurs: 



< 6 S < 15.79° 
15.79 <9 S < 18.34° 
18.34 <9 S < 39.26° 
39.26 <9 S < 43.71° 
43.71 < 0, < 44.67° 



Bistability of Rh S 2, m ,n and squares, e.g. (m,n) = (6,5). 

Rh s i,m,n stable, e.g. (m,n) = (ll,8). 

Everything unstable, e.g. (m,n) = (2,l). 

Rh S 2,m,n stable, e.g. (m,n) = (9,4). 

Super squares or anti-squares are stable, e.g. (m,n) = (19,8) 



44.67 < 9 S < 45.0° : Squares stable, e.g. (m,n) = (29,12). 



The results for 45° < 9 S < 90° are essentially the same with Rh s \ and Rh S 2 
interchanged. Recall that the size of the periodic box for which these results 



apply is given by v™ +n and that the aspect ratio of the rhombs (rectangles) 
is given by ^-^ for Rh slimi „ and by — for Rh s2 ,m,,n- in figure 4 we show an 
example of the bifurcation diagram close to onset for the case (m, n) = (19, 8). 
We know that either super squares or anti-squares are stable and have drawn 
the super square case. We have not calculated which of these two planforms 
is preferred at onset since this is determined at an 0(2(m + n) — 1) = 0(53) 
truncation of the bifurcation equations! 
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4-2 Case II: The degenerate case 73 = 5 — [3, K = +1 



In the degenerate case ■& = S— (3 the quadratic coefficient, e, is zero in equation 
(5) and both square and hexagonal lattices can give locally stable planforms. 
We have therefore evaluated the signs of the eigenvalues listed in both table 1 
and table 2 which are determined at a cubic truncation of the bifurcation 
equations. 

4-2.1 The stability of rolls, hexagons and squares 

Our results are presented in the form of a bifurcation set in the (/3, <5)-plane 
shown in figure 5. 

We first recap the relative stability results given in [1] for the subspaces Z4 = 
z 5 = z 6 = of equation (5) and z 3 = Z4 = of equation (4) before discussing 
the results of our extended analysis. First, in the subspace of the hexagonal 
bifurcation problem it is found: 

• Hexagons are stable in region 1. 

• Rolls are stable in regions 2,3 and 4. 

In contrast, for the subspace of the square lattice bifurcation problem: 

• Squares are stable in region 1, 2 and 3. 

• Rolls are stable only in region 4. 

Together the results for the two lattices suggest that in regions 2 and 3, for 
an unbounded domain, that hexagons are unstable to rolls but that rolls are 
unstable to squares. However, there is no formal way of directly studying if 
hexagons are unstable to squares. 

In our analysis of the finer lattices the results given above are modified. For 
a given hexagonal lattice we find hexagons and rhombs (Rh/j 3 mn ) are stable 
in region 1. Rhombs (Rh/j3 imn ) are stable in region 2 and rolls are stable 
only in regions 3 and 4. All other planforms are unstable. The position of 
the line dividing regions 2 and 3 is dependent on the choice of lattice. As the 
lattice angle approaches |, this line approaches the line a± — 0,2- The rhombs 
(rectangles) have aspect ratio depending on the lattice angle but lying between 
-tj and 1. 

For a given finer square lattice we find the stability results are unchanged 
for squares and rolls, however rhombs can be stable in regions 1 and 2. In 
particular, for the rhombs, we find that rhombs, Rh S 2,m,n are stable if < 
9 S < I and the rhombs, Rh s i )mj „ are stable if | < 9 S < |. These have aspect 
ratio dependent on 9 S , but again lying between -^ and 1. In addition, there are 
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regions of stable rhombs within region 1 which exist for | < 6 < sin 1 4^ « 
35.3° and for | > 9 > cos -1 4^ « 54.7° (aspect ratios between -^ and ^ ,~ ). 

In summary, although the regions of stability of the hexagons and squares are 
unchanged by the extended bifurcation analysis, the inclusion of the rhombic 
states allows for bistability of hexagons and rhombs of aspect ratio close to 
one. 



4-2.2 Unfolding the degenerate problem 

When e ^ the bifurcation equations for the hexagonal lattices contain a 
quadratic term and all planforms are necessarily locally unstable. However, if 
the quadratic term is sufficiently small compared with the cubic term, stable 
states may result through secondary bifurcations. For example, in the con- 
ventional analysis of hexagons using six critical wavevectors (figure 1(b)), for 
the degenerate bifurcation problem in regions 2,3 and 4 of figure 5, rolls and 
hexagons both bifurcate as pitchforks and rolls are stable. If a small quadratic 
term is added, then local to the trivial solution at /i c there are no stable 
bifurcating branches and the hexagons bifurcate transcritically creating two 
distinct branches, H + and H~. An example bifurcation diagram of such a sce- 
nario is shown in figure 6(a). In this figure it can be seen that as ft is increased, 
there is a hysteretic transition to hexagons from the trivial solution. On fur- 
ther increase in /z there is a second hysteretic transition, this time between 
hexagons and rolls. 

In our extended bifurcation analysis, again all planforms arise at pitchfork 
bifurcations in the degenerate case. When e ^ both super hexagons and 
hexagons now bifurcate as transcritical bifurcations and all states are locally 
unstable. Secondary bifurcations can, however, again stabilise some of the 
branches. If we consider j3 and 5 with values corresponding to region 2 of 
figure 5 and if n-cos 2 ^ ] is sufficiently small, then the bifurcation diagram 
of figure 6(a) is replaced by that shown in figure 6(b). The hexagons now 
undergo a hysteretic transition to rhombs rather than to rolls. Remember 
that the aspect ratio of the rhombs is dependent on the lattice: those shown 
in figure 6(b) are for the lattice (m,n) = (3,2) which results in rhombs of 
aspect ratio ^- ~ 0.87. 

The criterion that (1 e 2g ] is sufficiently small comes from requiring that the 
e corrections to the cubic coefficients are small enough to be neglected. If this 
is not the case, which is inevitable for fixed |e| <C 1 if the full range in 9h is 
considered, then a sequence of more complicated transitions can occur. We 
illustrate this by once again considering the case of Marangoni convection. 
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4-2.3 Unfolding the degenerate problem: the example of Marangoni convec- 
tion 

Although the Marangoni problem is nondegenerate, previous studies [13] have 
assumed that the quadratic term is sufficiently small and that it can legiti- 
mately be compared with the cubic terms. Specifically, the Marangoni prob- 
lem has f3 = —^7/8,5 = —3^/7/4,6 = -^- and although e does not appear 
to be small, the cubic coefficients themselves are relatively large. For exam- 
ple, ai = — ^ji an d a 2 = — ■j^- resulting in a saddle- node bifurcation on 
the hexagonal branch occurring at \x — —0.018. We have calculated all the 
eigenvalues that are determined at cubic order for the Marangoni case for the 
family of hexagonal lattices and found the corresponding bifurcation lines. The 
results are summarised in figure 7. The horizontal axis gives the lattice angle, 
9h- For the eigenvalues shown, this diagram is reflection symmetric about the 
line 9h — | and we have therefore only shown < 9h < f ■ For simplicity, 
we have only shown bifurcation lines which separate stable from unstable re- 
gions of the different planforms. For 9h approximately between 15° and 45° a 
transition from hexagons to rhombs occurs in a similar manner to that shown 
in figure 6(b). However, for 9h outside this range the e corrections to the cu- 
bic coefficients result in many further secondary transitions of significance. 
Five main shaded regions are shown indicating different combinations of sta- 
ble planforms, including regions of stable super hexagons with rhombs and/or 
hexagons for a cubic truncation of the bifurcation equations (note that the rel- 
ative stability of super hexagons and triangles is only determined at 0(2m — l) 
[8]). There is, in addition, a very narrow region to the left of the black region 
where super hexagons are the only stable planform. This region is too small 
to be readily discernible from figure 7, but is apparent in the bifurcation di- 
agram shown in figure 8. This bifurcation diagram was computed from the 
bifurcation equation (5) for the lattice (m,n) = (12, 7), i.e. 9h = 10.99°. Only 
branches which are stable for ji close to critical are drawn. In this particular 
case there is a region of bistability between hexagons and super hexagons. 
A hysteretic transition between hexagons and super hexagons can occur. An 
illustration of the super hexagon state for the case (m,n) = (12,7) is also 
shown. Although such complex patterns have not yet been seen in convection 
problems, recent experimental results of the Faraday crispation experiments 
do show such "superlattice patterns" [7]. 

Note that, if all lattices (all 9h) are considered, then hexagons are necessarily 
unstable (see < 9h < 8° in figure 7). However, in this case where we have 
fixed e 7^ 0, it is unclear whether this local problem is still valid for the 
amplitudes where secondary bifurcations arise. Nevertheless, it demonstrates 
some of the intriguing possibilities associated with this bifurcation problem. 
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5 Conclusions 



Standard low- dimensional bifurcation analyses of squares and hexagons give 
only a restricted stability analysis of these planforms. We have shown that, 
with one additional perturbation calculation, i.e. the calculation of C^, all 
the additional coefficients required to apply the extended stability analysis 
of [11] are determined. We have performed this calculation for the case of 
the long- wavelength convection equation (1) and analysed the results in two 
main cases. In Case I we found that extending the stability results significantly 
increased the known region of instability for squares in an unbounded domain. 
In particular, we find that squares are only stable at onset if e ~ 0. This is 
interesting given that it is already known that all solutions which are periodic 
on a hexagonal lattice are unstable at onset unless e = 0. For a given box 
with periodic boundary conditions we found that the predicted planform at 
onset was strongly dependent on the size of the box: some of the more exotic 
planforms such as super squares and anti-squares could be stable. In Case II 
we showed that regions of bistability of rhombs and hexagons exist. This gave 
a formal setting for studying the transition between hexagons and rhombs of 
aspect ratio close to 1 (although not squares). 
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(a) 
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Fig. 1. (a) Square lattice generated by four wavevectors on the critical circle oriented 
at 7£ to each other, (b) Hexagonal lattice generated by six wavevectors on the critical 



circle oriented at % to each other. 
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Fig. 2. (a) Square lattice generated by eight wavevectors on the critical circle. This 
case shows (m,n) = (2,1). (b) Hexagonal lattice generated by twelve wavevectors 
on the critical circle. This case shows (m,n) = (3,2). 
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Fig. 3. Bifurcation sets in case I for (m,n) = (2, 1). The lines indicate transitions 
in the number of negative eigenvalues, (a) Bifurcation set for rolls. The shaded 
region shows the region where rolls are preferred, (b) Bifurcation set for squares. 
The shaded region shows the region where squares are preferred. 
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Fig. 4. Bifurcation diagram for the square lattice (m, n) = (19, 8) for the Marangoni 
problem, f3 = —^C-,5 = — ^-. The stable branch is marked with a solid line and 
the unstable branches with dashed lines. An example of the super square planform 
is shown: this is spatially periodic, i.e. not a quasipattern. Only one period of the 
pattern is shown. 
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Fig. 5. Bifurcation set in the case e = 0. Only critical eigenvalues which divide stable 
from unstable planforms are shown. The different regions are discussed in the text. 
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Fig. 6. Unfoldings of the bifurcation diagrams relevant to region 2 in figure 5 for e 
sufficiently small, (a) The bifurcation diagram obtained by considering the subspace 
^4 = z§ = zq = of equation (5). (b) The new bifurcation diagram when the 
additional branches in our extended analysis is considered. Note that the rhombic 
(rectangular) pattern has an aspect ratio close to f. 
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